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Abstract 

Based on a technique of Barvinok (H |5l [6l and Barvinok and Soberon (8l |3 we 
identify a class of edge-coloring models whose partition functions do not evaluate to 
zero on bounded degree graphs. Subsequently we give a quasi-pol 5 momial time ap¬ 
proximation scheme for computing these partition functions. As another application 
we show that the normalised partition functions of these models are continuous with 
respect to the Benjamini-Schramm topology on bounded degree graphs. We moreover 
give quasi-pol 5 momial time approximation schemes for evaluating a large class of graph 
polynomials, including the Tutte pol 5 momial, on bounded degree graphs. 

Keywords: Benjamini-Schramm convergence, approximation algorithms, partition functions, 
sparse graph limits. 


1 Introduction 

In a series of papers Barvinok IK [5l 0 and Barvinok and Soberon (9| developed a 
technique that yields quasi-pol 5 momial time approximation schemes for several hard 
counting problems such as computing the permanent of a matrix and computing the 
number of (edge-colored) graph homomorphisms into a fixed graph. Their technique 
is essentially based on two things. First the counting problem is cast as the evalua¬ 
tion of some pol 5 momial and they identify a region where this pol 5 momial is nonzero. 
Secondly, they prove that in this zero-free region the logarithm of the pol 5 momial is 
well-approximated by a low order Taylor approximation, which they prove to be com¬ 
putable in quasi-polynomial time. 

In this paper we will apply this method to partition functions of edge-coloring 
models and more generally to Holant counting problems and contractions of tensor 
networks. That is, fixing a degree bound A, we will identify edge-coloring models 
for which the partition function does not vanish on graphs with maximum degree at 
most A and consequently show how to approximately compute the partition function 
in quasi-poIynomial time. We moreover show that the normalised logarithm of the 
partition function is continuous with respect to the Benjamini-Schramm topology for 
these models, and we give give quasi-poIynomial time algorithms for approximately 
evaluating a large class of graph polynomials, including the Tutte polynomial and the 
matching polynomial, on bounded degree graphs. 

Below we will introduce and describe these concepts in a bit more detail and state 
our results. 

* University of Amsterdam. Email: guusregts@ginail.com. 
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1.1 Partition functions of edge-coloring models 

Edge-coloring models originate in statistical physics and their partition functions have 
been introduced to the graph theory community by de la Harpe and Jones 11201 (where 
they are called vertex models). We call any map h : t C a k-color edge-coloring model. 

(In this paper we will fix A: > 2 and just speak about edge-coloring models). For a graph 
G = {V,E), the partition function of h is defined by 

P(G){h):= ^ ( 1 ) 

<p:E^[k]veV 

where S{'v) denotes the multiset of edges incident with the vertex v and where (p{S{v)) 
denotes the multiset of colors that the vertex v 'sees', which we identify with its inci¬ 
dence vector in so that we can apply h to it. Let us give some examples. 

Example 1. • Let h : —)• C be defined by h{cii,ci2) = 1 if aj < 1 and zero 

otherwise for { oli , 0.2) G N^. Then p{G){h) is equal to the number of matchings 
of the graph G. This can be easily seen as follows: the edges receiving color 1 
correspond to a matching if and only if the coloring gives a contribution to the 
sum (1T]|. 

• Fix d G N and let h : —)• C be defined by h{a.i,0.2) = 1 if ai = d and zero 

otherwise for (ai,a 2 ) G N^. Then p{G){h) is equal to the number of d-regular 
subgraphs of G on the same vertex set. 

See |[37 l and below for more examples. 

Interchanging the role of edges and vertices in |[T]| one obtains partition functions 
of vertex-coloring models. Let a G C” such that a, 7 ^ 0 for all i and let B G be 

a symmetric matrix. We call the pair (a, B) a vertex-coloring model (In Il20i this is called 
a spin model provided a is the all ones vector). For a graph G = {V,E) the partition 
function of (fl, B) is defined by 

p(G)(fl, B). n n ^ 

(p:V^ln]veV uveE 

If a is the all-ones vector and B is the adjacency matrix of some graph H, then p(G) {a, B) 
is equal to the number of graph homomorphisms from G to H. 

Partition functions of edge-coloring models include partition function of vertex¬ 
coloring models, as has been shown by Szegedy l37ll38l , cf. Lemmaj^in Section lZTI But 
they form a much bigger class of graph invariants, cf. II2TI . For example, the number of 
matchings is not equal to the partition function of any vertex-coloring model. Partition 
functions of edge-coloring models have been characterised in l37l and ll2TI based on 
invariant theory of the orthogonal group and the Nullstellensatz. 

In many cases it is NP-hard to determine whether or not the partition function of an 
edge- or vertex-coloring model vanishes on general graphs. However, when restricted to 
bounded degree graphs things may become easier. In particular, Barvinok and Soberon 
[18| proved that when a is the all ones vector and B is sufficiently close to the all ones 
matrix, then 


p(G)(fl, B) 7 ^ 0 for all graphs of maximum degree at most A. (3) 

Based on the technique of Barvinok lH |5l | 6 | and Barvinok and Soberon ( 8 ) |9l we will 
prove in Section a similar result for partition fimctions of edge-coloring models: 

Theorem 1 . Let A G N then there exists a constant f > 0.35 such that for any k-color edge¬ 
coloring model h that satisfies \h{a.) — 1| < /3/(A -|-1) for all a. G with |a| < A we have 
for any graph with maximum degree at most A, p(G) {h) 0. 


2 


In fact, we will state and prove in Section a slightly more general version of The¬ 
orem [1] that also applies to tensor networks. The definition of tensor networks will be 
given in Section |2l 

Remark 1. While every partition function of a vertex-coloring model is equal to the 
partition function of some edge-coloring model, the exact relation between Theorem [T] 
and @ is presently not clear. It is possible though to derive a qualitative version of 
(|3} from Theorem [T] (even valid for vertex vertex-coloring models with vertex weights 
close to 1). We will however not do this here. Instead, in Section |3] we will present an 
example of a vertex-coloring model to which the result of Barvinok and Soberon does 
not apply, while from our results it follows that its partition function does not evaluate 
to zero on bounded degree graphs, cf. Example lU 


l. 2 Computing partition functions 

Exact computation Computing partition functions is in many cases #P-hard. There is 
a line of work by Dyer and Greenhill 1221 , Bulatov and Grohe 1131 , Gai, Ghen and Lu 

m, where dichotomy theorems have been proved for the complexity of computing 
partition functions of vertex-coloring models. Roughly they proved that if the model 
has a very special structure, its partition function can be computed in polynomial time 
and it is #P-hard to compute otherwise. In a sequence of papers Gai Lu and Xia IIT7l[T8ll , 

Gai, Huang and Li 11611 , and Gai, Guo and Williams IITSi have proved similar results for 
partition functions of edge-coloring models with 2 colors^ 

Approximate computation Since computing partition functions is in many case #P- 
hard, people have started looking at approximation algorithms. Let us consider the 
case of counting weighted independent sets. Lor a graph G and A > 0, denote by Pa{G) 
the sum over all independent sets of G, where an independent set I is counted with 
weight aI^I. Note that p\{G) may be viewed as the partition function oi a = (A, 1) 

and B = ^ ^ | ^. A surprising threshold phenomenon has been found for the com¬ 
plexity of approximating pa{G): for A < Ac(A) := (A — (A — 2)^^ there exists a 

fully-polynomial time approximation scheme (EPTAS) for graphs of maximum degree 
at most A, discovered by Weitz 1391 ; for A > A^ there does not exist a fully-pol 5 momial 
time randomised approximation scheme (EPRAS) for A-regular graphs (with A > 3), 
unless NP=RP, as has been proved by Sly and Sim 1361 . The threshold Ac is the critical 
point for the (infinite) A-regular tree to have a unique Gibbs measure. This phenomenon 
suggests an intimate cormection between phase transitions in statistical physics and 
computational complexity. The ideas of Weitz have been applied and extended to other 
types of counting problems and partition functions; see e.g. l3l 12411311 l35l . Other in¬ 
teresting approximation algorithms and hardness results have been found for partition 
functions of vertex-coloring models an planar graphs l26l . 

Our contribution In Section |4] of this paper we apply Theorem [1] to obtain a quasi¬ 
polynomial time approximation scheme (QPTAS) for the computation of partition func¬ 
tions of a large class of edge-coloring models on bounded degree graphs, based on 
a technique developed by Barvinok |4l]. This technique was used by Barvinok and 
Soberon |3 to obtain a QPTAS for the computation of partition function of vertex¬ 
coloring models on bounded degree graphs (with a equal to the all ones vector). 

The main idea is to rewrite the partition function as the evaluation of a univariate 

^In fact, in these papers they consider Holant counting problems, but they include partition functions of edge¬ 
coloring models. These Holant counting problems may be seen as tensor network contractions and we will 
define these in Section lA2l 
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polynomial q at 1 and then use a low order Taylor approximation to In q. The approxi¬ 
mation turns out to be good when the roots of q have absolute value at least M for some 
fixed constant M > 1. Theorem [T] is used to guarantee the existence of this constant M. 
From the perspective of the Lee-Yang approach to phase transition in statistical physics 
Il28l , one may view zero-free regions of the partition function as regions far away from a 
phase transition. So this approach also indicates a connection between statistical physics 
and computational complexity. 

We need one definition to state our result. 

Definition. Let e > 0 and a e C. We call b e C a multiplicative e-approximation to a if 
and if the angle between a and b is at most e. We call b e C an additive 
e-approximation to a if |fl — b| < e. 

We prove the following result in Section ID 

Theorem 2. Let A e N. Then there exists a constant <5 > 0 (we may take 5 = 0.35/(A -|- 1)1 
such that for any e > 0, any edge-coloring model h that satisfies \h{a) — 1| < 5 for each a e 
and any graph G = {V,E) of maximum degree at most A, we can compute a multiplicative 
e-approximation to p{G){h) in time |and an additive e\V\-approximation to 
In |p(G)(b)| in time | 

Theorem|2]also applies to contractions of tensor networks and Holant counting prob¬ 
lems. See Section 3] for the exact details. 

1.3 Evaluating graph polynomials 

It turns out that the technique of Barvinok can also be applied to graph polynomials 
of exponential type. This is a class of graph pol 5 momials introduced by Csikvari and 
Frenkel 1T91 . which includes many important pol 5 momials such as the Tutte polynomial, 
the matching polynomial and many more. We will define these polynomials in Section 
14.31 where we obtain a QPTAS for computing certain evaluations of these pol 5 momials 
on bounded degree graphs, cf. Theorem IT^ Let us for now specialise these results to 
the Tutte pol 5 momial. 

For a graph G = (V, E) the random cluster model formulation of the Tutte polynomial 
Z{G){q, v), is defined as follows: 

Z{G){q,v) := (4) 

ACE 

where k{A) denotes the number of connected components of the graph {V,A). We will 
consider the Tutte polynomial for fixed z; as a pol 5 momial in q. 

Theorem 3. Let A e N and let vq E C. Then there exists a constant C = C(A, i^o) such that 
for any e > 0, any graph G = {V,E) of maximum degree at most A and any q E C with |(j| > 
C, we can compute a multiplicative e-approximation to Z{G){q,VQ) in time and 

an additive e\V\-approximation to In \Z{G){q,Vo)\ in time 

In l25l l Goldberg and Jerrum showed that approximating Z(^G){q,v) for general 
graphs and for q > 2 and z? > 0 is as hard as counting independent sets in bipar¬ 
tite graphs (#BIS). The present results may be an indication that this is not the case for 
bormded degree graphs. See 1311 for related algorithmic results for the Potts model and 
1231 for hardness results for the Potts model on bounded degree graphs. 

1.4 Sparse graph limits 

In 2001, Benjamini and Schramm flOl defined a convergence notion for finite bounded 
degree planar graphs, which can also be extended to general bormded degree graphs 
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and is often referred to as Beniamini-Schramm convergence, local convergence or weak con¬ 
vergence. This convergence notion has many equivalent definitions (see the book by 
Lovasz Il2^ Proposition 5.6]) and we choose one which is most suitable for our pur¬ 
poses. Fix a degree bound A and let (G,,) be a sequence of simple graphs, each of 
maximum degree at most A. We call the sequence (G,,) locally convergent if \V{G„)\ 
goes to infinity, as n —> oo and if for each cormected simple graph H the quantity 


md(H, G„) 

\y{Gn)\ 


is a convergent sequence of real numbers, where in ||5ll ind(H, Gn) denotes the number 
of induced subgraphs of G„ that are isomorphic to H. As an example let C„ be the cycle 
of length n. Then (C^) is locally convergent. See |29| for possible representations of 
limit objects. 

Let be the collection of simple graphs of maximum degree at most A. Call a 
graph parameter / : > C estimable if for any locally convergent sequence of graphs 

(Gn) C Ga, (/(Gn)) is a convergent sequence of numbers. In Section |5| we look at the 
estimability of the normalised partition function n{G){h) of an edge-coloring model h, 
which is defined as 


n{G){h): = 


ln|p(G)(/z)| 

I^G)! 


( 6 ) 


where G is a graph. This normalisation is motivated by statistical physics, where this 
quantity is known as the entropy per vertex. Note from the perspective of estimability it 
does not matter to take the absolute value of p{G){h) before applying the logarithm. For 
if we fix a some branch of the (complex) logarithm, then the imaginary part is bounded, 
and dividing it by |y(Gn)| will make it disappear in the limit. 

Borgs, Chayes, Kahn and Lovasz KTTI showed that if one in ® replaces p{G){h) 
by p{G){a,B) for appropriate (fl,B) this parameter is estimable; see also l30ll . In 119^ 1 
Csikvari and Frenkel showed that if one replaces p{G){h) by q{G){z) in where q 
is a graph pol 5 momiaI of bounded exponential type and |z| is large enough, then this 
parameter is estimable. This extends work of Abert and Hubai |T|| who showed that this 
holds for the chromatic polynomial. 

Utilising Theorem [Tj and building on a result of Csikvari and Frenkel lfT9l , we will 
prove that the normalised partition function of an edge-coloring model that is suffi¬ 
ciently close to the all ones edge-coloring model is also estimable. 


Theorem 4. Let A e N. Then there exists a constant f> > 0.35 such that for any edge-coloring 
model h that satisfies \h{a) — 1| < f/{A + 1) for each a e N/ the normalised partition 
function n{-){h) : Ga —t IR is estimable. 


1.5 Organisation of the paper 

In the next section we collect some preliminaries and set up some notation. In particular, 
we introduce tensor networks. In Section we will prove a slight generalisation of 
Theorem [1] valid for tensor networks. In Section |4| we will discuss the technique of 
Barvinok |4j mentioned earlier and prove Theorems |2| and |3j In Section|5|we will give a 
proof of TheoremlH Finally in Section[^we will conclude with a few remarks. Sections 
[Hand 13 may be read independently of each other. 
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2 Preliminaries 


2.1 Notation and definitions 

We will setup some basic notation and give some definitions here that are used through¬ 
out the paper. 

By Q we denote the collection of all graphs, allowing multiple edges and loops. 

Recall that by we denote the collection of simple graphs of maximum degree at 
most A. Let S be a ring. An S-valued graph invariant or S-valued graph parameter is a map 
f : Q ^ S that maps isomorphic graphs to the same element of S. We often omit the 
reference to S. When S = C[z] (the polynomial ring in the variable z) we call / a graph 
polynomial. A graph invariant f is called multiplicative is /(G U H) = f{G)f{H) where 
GU H denotes the disjoint union of two graphs G, H. Clearly, partition functions of 
edge-coloring models are multiplicative. 

Whenever we talk about an edge-coloring model we will always assume it has k 
colors, where k is some fixed natural number. In particular, in algorithmic statements 
we will consider /c as a constant and moreover assume that it is at least 2, for else 
computing the partition function is trivial. 

Let U be any kx n matrix and let a G C". Let mi, ..., be the columns of U. Define 
the fc-color edge-coloring model h^ jj by 

n k 

,=1 ;=1 

for a G In case a is the all ones vector we just write hu instead of ha jj. Szegedy 
Il37l l38l showed the following 

Lemma 5. Let B be an n x n complex symmetric matrix and let a G (C \ {0})”. Let U be 
any (complex) matrix such that U^U = B. Then for any graph G we have p{G){a,B) = 

p{G){Ku)- 

2.2 Tensor networks 

As mentioned in the introduction. Theorem [T] not only applies to edge-coloring models 
but also to tensor networks. We will introduce tensor networks here. 

First we need some further definitions. For a G set |a| := Ef=i By we 

denote the set {a G : |a| = d}. We identify the set of functions {/z : t C} with 

and denote it by S^. Let G = (y, E) be graph. We define 

= n 

v^V 

Let h = (E^jugv G Sg- We call the pair (G,E) a tensor network. The contraction of the 
tensor network {G,h) is defined as follows: 

P(G)(E):= ^ (8) 

<p-.E^[k]veV 

This is almost the same as the partition function of an edge-coloring model. The main 
difference is that in the tensor network we can select a function for each vertex, while 
for an edge-coloring model the function is only allowed to depend on the degree of 

^ While Szegedy only proved this for real matrices B, it is easily seen that the proof remains valid for complex 
matrices, cf. l33l Lemma 7.1]. 
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the vertex. It will be convenient to view an edge-coloring model h as an infinite vector 

The reason why (G,h) is called a tensor network comes from fhe facf fhat we may 
view as a symmefric tensor in for each v ^ V. The confraction of the 

tensor network then corresponds to contracting a tensor obtained from fhese h^’s along 
the edges of G with respect to a bilinear form. We refer to Section 3.4 of fhe aufhor's 
PhD fhesis 1331 for more defails on fhis perspecfive. 

Confractions of tensor networks have been used to simulate quantum computing by 
Markov and Shi in 132| , where they gave a polynomial time algorithm for confracting 
tensor networks on bounded degree graphs of bounded freewidth. In {2] Arad and 
Landau use quantum algorithms to compute additive approximations to tensor net¬ 
work contractions. We refer to both l32l and [2] and the references in there for more 
defails concerning the connection of fensor networks and quantum information theory. 
Contractions of fensor networks can also be used to model Holant counting problems 
(see CZllIIllTHl). 

2.3 The orthogonal group 

There is a natural action of fhe complex orfhogonal group, Oj-(C) (the group of k x k 
matrices g such that g'^g = I where I is the identity matrix) on tensor networks (and 
hence on edge-coloring models) that leaves the tensor network contraction invariant. 
To define fhis acfion in a convenient way, we need to realize that {h : > C} is 

in one-to-one correspondence with the space of linear map from C[xi,.. .,xj;] to C. 
Each such map is determined by its action on monomials and the monomials are in 

one-to-one correspondence with This we may view any h 6 C'’^^ as a linear map 
h : C[xi,..., Xjt] —> C that is zero on all monomials of fofal degree bigger than d. Now 
Oi-(C) acts on C[xi,...,xj-] as follows: an element g e 0/f(C) sends a polynomial p to 
the rmique pol 5 momial gp that on input v ^ takes the value p(g'^v). Then Oi((C) 
also naturally acts on the space of linear maps C[xi,.. .,Xi^] —> C via gf(p) := fig^p) 
for g e 0/c(C) p e C[xi,..., Xj^], V e and any linear map / : C[xi,...,xj-] —>■ C. 

For a fensor network {G,h) with h = e Sq and g e Ojt(C) we define {G,gh) 

fo be the tensor network with gh := {gh^)vGV G Sq- Then 

for any g e Oi.(C), graph G and h e Sc- p{G){gh) = p{G){h). (9) 

For edge-coloring models this is proved in {21] and the same proof also works for 
confractions of fensor networks. 


3 Absence of roots for tensor network contractions 

In this section we will find zero-free regions for fensor nefwork confractions. In partic¬ 
ular, Theorem [1] will be proved here. We begin with some notation, necessary to state 
our main theorem. Define for d e N, <5 e (0,1) and tj > 0, the set 

> rj, for all a, a' e 


CSg- 


Sd{8,rj) :={/!eC^S | \h^-h^,\ <5 aod\h^\ 
and for a graph G = {V,E), fhe sef 


We can now sfate our main fheorem. 
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Theorem 6 . Let G = (V,E) he a graph and let rj e (0, oo). For any 0 e {0,2n/3) fix 
f> < rj9cos{9/2) and let 5 := min{rj,. Then for each h E Sc{h,r]), \p{G){h)\ > 
(cos(0/2);/)l'^lfcl^L In particular, p(G)(h) 7 ^ 0. 

We will prove Theorem [ 6 ] in the next section, but we will first state some conse¬ 
quence. 

Obviously, the goal is to choose 9 with f as large as possible. For this we need 
to solve maxgg(Q 27 r/ 3 ) ^cos(0/2). This maximum is attained by the solution 9* to the 
equation 2/9 = tan(0/2), which is approximately equal to 9* w 1.72067. Let x* = 
9* cos(9*/2). Then x* w 1.12219. Define for d e N 


r(d) 


X 


* 


1 + 


X* 

23 


Then 0.71885 w ^*(1) < f*{2) < f* (3) <■■■< x* ^ 1.12219. We have the following 
corollary, which implies Theorem [1] 

Corollary 6a. Let G = {V,E) be a graph. Fix nonzero a^ E C for each v E V. Write 
f>* := f*{A{G) + 1). Let h E Sq such that \av — h^{ix)\ < 2 {a{g)+i) all v E V and 

“ e ^degivy 

|p(G)(fz)| > n l«-l ■ (cos(r/ 2 )(l - 2 a(g)+ 2 ^)"''^'"' 

In particidar, p{G){h) 7 ^ 0. 

Pt* 

Proof. Let us first consider the case where fly = 1 for all v. Lef t] := 1 — 2 a{g)+ 2 ' ^ 

o* 

simple computation shows that f>* = rjx*. Clearly, t] > f>* /2. So if we let 5 = 
then each E S{S, rj) and hence ifTOt . with fly = 1 for each v eV, follows from Theorem 

m 

In the general case let g^{oc) = h^{a) /fly for each a and v E V. Then for each v and 
|1 “ ^ 2 (a(g)+i) ' flz; = 1 for each v eV, holds for p{G){g). Since 

p{G){h) = riygvfli; ■ p{G){g), lO now follows. □ 


We obtain another corollary if we take the action of the orthogonal group into ac¬ 
count. Let X e C^ and consider x as a /c x 1 matrix. Let hx denote the edge-coloring 
model as defined in lO (so n = 1, u^ = x and aj = 1 in this case). 

Corollary 6 b. Let G = {V,E)be a graph. Then for each x E such that (x, x) = xj 7 ^ 

0 and nonzero ay E C for each v E V, there exists a constant d > 0 (only depending on x and 
A(G)) such that ifh E Sg satisfies 


\h-ix) 


avhxi‘>i)\ < 


fly||(X,X)/fc|‘i‘=gW/2^ 
2(A(G) + 1) 


for ally eV and tx E then p(G)(h) 7 ^ 0. 

Proof. Since (x, x) 7 ^ 0 there exist g E Ojt(C) such that gx = Al, a scalar multiple of the 
all-ones vector, with A = ((x,x)/A:)^'^^. Then for any a e we have fhat (a) = A'^. 
By continuity, if h^ is close to ayhx it follows thaf gh°{oi) is close to flyA*^ for each a E 
The resulf now follows from Corollary and the fact that tensor network contraction 
is invariant under the action of the orthogonal group. □ 
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Corollary |6 ]d is roughly saying that if the tensors are close enough to the rank one 
tensor hx, then the tensor network contraction is not equal to zero. 

We will now give an example illustrating the power and weakness of our results. 

Example 2. Consider for A, ^ € C the matrix ^ ^ ^ ^ ~ (^' ^)- Since 

the matrix B contains a zero, the result of Barvinok and Soberon [|8| does not apply to 

it. However, if we let U = 1/ \/2 f then U^U = B. Consider the edge- 

coloring model hjj as defined in |(7ll. Then, by Lemma|^ we have p(G) (fl, B) = p(G) (hjj) 
for each graph G. By Corollary [6^, for fixed nonzero p, fixed A > 0 and small enough 
A we have that p(G)(lz(j) t^O for all graphs with maximum degree at most A. 

For p = 1 we can view p(/i(j)(G) as a polynomial in A, which for A-regular graphs 
coincides with the independence polynomial of G evaluated at A'^. Scott and Sokal f34 l 

proved that, if |A^| < —, fhen p{hu)[G) ^ 0 on any A-regular graph G. They 

moreover showed thaf this bound is tight. As A —> oo this boimd behaves roughly like 
while the result we obtain is of the order O(^), showing that our results are not 
optimal. 

3.1 A proof of Theorem [6] 

Our proof of Theorem is based on a method of Barvinok aEiiii and Barvinok and 
Soberon |8l|3. Following Barvinok |0 we can say that although the method of proof is 
similar to the methods considered in HI 018] |3 O it does require some effort and new 
ideas. 

We will first state and prove several lemmas. The first lemma about angles between 
vectors is due to Boris Bukh see |5j . We refer fo ISj for a proof. 

Lemma 7. Let i9 E [0,27r/3). Let be nonzero vectors in such that for each pair 

1.1 the angle between x, and Xj does not exceed &. Then || x,|| > cos(d/2) ||^!|l- 

We will consider elements of C as vectors in IR^ so that we can speak about the angle 
between two complex numbers. This is how Lemma [7| will be applied. 

It will be convenient to introduce some notation. Let F, F' be finite sets with F C F'. 
Let ^ : F —t [fc] and let xp : F' ^ [k]. We say that ip extends cp if xplp = cp. When G = {V, E) 
is a graph, u E V, F C E and a E we say that a is u-compatible with (p if there 

exists a map cp' : F LI 3(u) —)• [k] that extends cp such that (p'{S{u)) = a. We denote this 
by f ci. 

It will also be convenient to consider the coordinate ring of C'^rf, which is the poly¬ 
nomial ring := C[xa | a € IN^]- For a graph G = {V, E) we define 

Rg ■= n ^degivy 

17€ V 

For F C E and <p : F ^ [k] define 

P^(G):= E (11) 

Tl>:E^[k]veV 

f\p=<p 

So p^(G) is a polynomial and we can evaluafe this at h E Sq. In particular, if F = 0, 
then p^{G){h) just coincides with the ordinary tensor network contraction p{G){h). 

In the two lemmas below we assume that we have fixed some graph G = (V, E) and 
p, <5 > 0. Moreover, for u E V we let (5 (m) C E be the set of edges incident with u and 
N(w) C V the set of neighbours of u. 
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Lemma 8 . Let t > 0, let F C E, let (p : F ^ [fc] and let u E V. If for all h e Sc{S,r]) and 
for each xp : FU S{u) —> [k] extending cp we have 7 ^ 0 and moreover for each 

V e N{u)U {u}: 




aelN^ , - 

deg(c) 

ip^vOi 


dxl 




( 12 ) 


then the angle of p^^^^‘\G){h) and 
(p is at most 


p^^^^‘\G){h) for any xp,xp' : F\jS{u) —> [k] extending 

Hm±2i. , 13 ) 

TTj 


Proof. Let us write F' = FG5{u) and fix \p : F' —r [k]. Since p^'(G) 7 ^ 0 on Sq{S,p) 
we can consider a branch of ln(p^^(G)) on Sq{5,p). For any v e N(m) U {u} and 

^ ^deg(i;) 




ln(pf(G)) 


(^) 

P^'(G) 


(14) 


Note that (IT4ll equals zero if a is not i;-compatible with xp. Then for each h e Sg{S,p) 
we obtain by (IT2ll that 


aelN^ , , 

deg(D) 


^Mp?(g)(p)) 


< 


Tmin{|P„| I h e Sg{S,p),x e N^eg(i;)} 


1 

Tf] 


< —■ (15) 


Next, let xp' : F' ^ [k] be another map that extends xp and define an element g = 
e Sg{S,p) as follows. If v i N{u) U {u}, then gl = hi for all a e N^eg(z;)- ^ 
V e N{u), lef E{uv) denote the set of edges cormecting u and v (G mighf have multiple 
edges). Consider f = xp(E(uv)) and f' = xp'[E[uv)) as elements of N*^. Then for each 
DC e set a' := a — /3 + f. Then define gl := hi, if a' > 0 and gl = hi otherwise. 

Finally, set gl = /j“ if a 7 ^ xp{5{u)) and gl = otherwise, where f = xp'{S{u)). Note 
that p^,{G){h) = p^{G){g). This implies by the fundamental theorem of calculus thaf 
|51n(p^^(G)(P)) — 51n(p^)(G)(/i))| is bounded by 


max 

deg(n) 




\g:-K\< 


S{A{G) + l) ^ 

xp 


where 5/ denotes the imaginary part of a function / and where the inequality is due to 
(HI and the fact that \hl — hl,\ < S for any a, a' G ^deg(y)' I'l'l® proves the lemma. □ 

Lemma 9. Let 6 G [0,27r/3). Let w G V and let F C E with S{u) C F. Let cp : F ^ [k]. 
Suppose that for each v G N{u) U {u}, all h G SG{d,ij) and all maps xp,xp' : FU<5(x) —> 
[P] extending cp, we have p^^^^\G){h) 7 ^ 0 and that the angle between p^^^^\G){h) and 

p^^^^\G){h) is at most 9. Then for each v G N{u) U {u} and all h G Sg{S,p) we have 


|p^(G)(P)| > cos(F/2) \K\ ^pl{G){h) 


aeN^ , , 

deg(ii) 

<p~vOi 


(16) 


10 














Proof. First note that ifl^ is trivially true for v = u (since we could replace cos(0/2) 
by 1). Now let v e N{u) and let a e ^deg(D)‘ u-compatible with (p, then 

^P^{G) = 0; if a is u-compatible with cp, then 


9 p ^ 1 


shd(G) = E pr''"’(G). 


Since 


Tp:FUS{v)—^[k] 

ip\p=<p,xp{S{v))=a 


FUS(v)f 


we obtain that for any h e SQ{S,t]), 


pliG) = 

= E 



ip:FUS{v) 

^[k] 


i’\F=4 


cos(0/2) 

E 



tp:F\JS{v)^[k] 

ip\F=<p 

cos(0/2) 

E 



deg(B) 

4’ 

4’\f 



cos(0/2) 

E 

K\ 


FUS(v), 


E pr'“N)(/o 




f:FUi}(v)->‘[k] 


i,4icm 


deg(B) 


(17) 


(18) 


where the first inequality in iflSll is due to Lemma [71 since by assumption for each ip, ip' 
and any h e Sc{5,rj) we have that the angles between and 

differ by no more than 6, the second inequality is obvious and the equality in (ITSl is 
due to (ITTll . This proves the lemma. □ 

We can now give a proof of Theorem [6| 

Proof of Theorem\^ Fix a graph G = {V, E) and h e Sg{S, rj). If we plug in t = cos(0/2) 
in |[T3] |, then by assumption we have that dTSl l is at most 0. We will first show that the 
following three statements hold: 


(i) 

(ii) 


(hi) 


For each F C E, for each (p : F ^ [k] and for all E e SQ{5,rj) we have 

vli.G){h)fO. 

For each F C E and each m 6 V let : F —> [k] and write F' = FU S{u). 
Then for each xp, ip' : F' —>■ [k] extending cp and each h e SciS,?]) 
the angle between p^{G){h) and p^,{G){h) is at most 0. 

For each F C E and u Fz V such that S{u) C F let ^ : F —>■ [E]. 

Then for all h e Sg(<^, p) and each v e N{u) U {u} we have 


\pl{G){h)\> cos{e/2)E,^^^ 

^ deg(o) 

<p~vOi 


\K\ ^y^GUh) 


(19) 


Note that (fl^ (i) for F = 0 already implies that p{G){h) 0 for all ki G Sg{S,p). 

We will now prove that ifl^ (i), (ii) and (hi) hold for all sets F C E by induction 
on |E \ F|. If F = E, then (i) holds since p^{G){h) is just the product of |y| nonzero 
numbers. Clearly, (ii) also holds. Moreover, (hi) also holds since for each v, 


E K^pl{G){h) = p^{G){h) 

aelN^ , , “ 

deg(i;) 
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and the sum in this equation consists of only one term. 

Let now F be a subset of E of size strictly smaller than E. Let u e V. If C F, 
then (ii) is clearly true. So we may assume that S{u)<^F. Let (p : F ^ [k] and let 
: F U S{u) —> [E] be maps that extend (p. Then by induction from (i) and (iii) 
(since |FU(5 (m)| > |F|) we know that the conditions of Lemma 0 are satisfied with 
T = cos (0/2). This implies that (ii) holds for the set F. 

To prove that (i) holds for F, fix v such that S{v)^F. Then 

P^(G)(^)= E 

ip:FljS(v)—>-[k] 

1’\p=‘P 

As by induction the are nonzero and by (ii) we know that their angles 

differ by at most 9, Lemma [7] implies that p^(G)(h) is nonzero for any h € SQ{S,tj). 

Finally, take u such that S{u) C F. If no such u exists there is nothing to prove 
for (iii). Fix v G N(m)U{m}. Since the sum in ifl^ has only one contributing factor if 
S{v) C F, we may assume that 5{v)’^F. By induction, from (ii) we know that the angle 
between p^^^^^\G){h) and p^^^^\G){h) for any ip,ip' : FU<5(i;) —> [k] extending (p is 
at most 0,as |FU<5(f)| > |F|. So Lemma implies that (iii) holds for F. So by induction 
we conclude that ([19) holds for all sets F. 

To finish the proof write V = {vi,... ,v„} and let Fq = 0. For i > 0 let F; : = 
Fi_iLI S{vi). We claim that for each f = 0,... ,n, (p : Fi —)• [k] and E G Sc{S,p) we have 

|pJ(G)(E)|>fcl^\^>l<(cos(0/2))”-‘. (20) 

Indeed, for i = n this is clearly true. Next, pick any i < n and (p : Fj ^ [k]. Then 


\p^^{G){h)\= E P^HG)(P) >cos(0/2) E IP^nG)(E)| 
i’lf. =‘P '/’Ip. =<P 

> cos(0/2)/cl^'+i\hlfcT\h+il cos(0/2)"-('+i)p”, 

where the first inequality follows from (IT^ (ii), and the second by induction. This shows 
(120) and finishes the proof. □ 


4 Approximation algorithms 

In this section we will consider algorithms for computing tensor network contractions 
and evaluations of graph pol 5 momials. Both algorithms are based on a method of Barvi¬ 
nok [4J (which has also been used by Barvinok j|5l|6) and Barvinok and Soberon (8[|2), 
which we describe first. 

4.1 Approximating evaluations of polynomials 

Let p G C[z] be a polynomial of degree d and suppose that p(z) ^ 0 for all z in a closed 
disk of radius M. In this section we will describe the method of Barvinok |Jj to find 
a multiplicative approximation to p(f) for f G C, which is only based on knowing the 
values of ^p|t=o- 

Define the univariate function / by 

/(z):=lnp(z) (21) 
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for 2 G C, where we fix a (continuous) branch of logarithm by fixing the principal 
branch of the logarithm for p(0). 

To approximate p at f G C, we will find an additive approximation to / at f us¬ 
ing the the Taylor expansion around z = 0. This can then be transformed to give a 
multiplicative approximation to p. 

Let 


Tn(f)(t) :=/(0)+ £ 

m—1 


— — ffz) 
m\ dz^ 


z=0 


( 22 ) 


Note that (l22ll only depends on the values of ^p\t=o for m = 0,... ,n. (Where we agree 
that the 0*^ derivative is just the function itself.) To see this note that f'{z) = p'{z) /p(z), 
that is, p\z) = p{z)f{z). So for m> 1 we have 




t=o 



(23) 


This implies that if we can compute the values of ^p(z)|z=o for j = 0,.. .,n, then 
(l23] | provides a nondegenerate (as p(0) 7 ^ 0) triangular system of equations to compute 
^/(z)|z=o for i = 1,...,n, which can be done in time To summarise: 

Lemma 10. If the values ■^p\z=o are given for m = 0,... ,n, then ^f\z=o can be computed 
in time for each m = 0,... ,n. 

The quality of the approximation (l22ll depends on the complex roots of p. 

Lemma 11. Let p be a polynomial of degree d such that p(0) 7 ^ 0/or all z such that |z| < M 
for some M > 0 and let t & C such that |f| < M. Then for n = 0(lnl/e) ive have 

||T„(/)(f)|-ln|p(f)|| <d 6 , (24) 

and for n = 0{\nd/e) we have 


< 


eMfm) 


P(f) 


< e‘^ and the angle between p(f) and is at most e. (25) 


Proof. Let q := \ t\/M. We will first show 

l/(0-T„(/)WI< („//(;'_^) . (26) 

The proof of ll26ll is quife similar fo fhe proof of Lemma 1.2 from | 8 j, but for complefe- 
ness we will give if here. 

Since p(0) 7 ^ 0 we may write 

p(z) = p(0)n(^l-|), (27) 

where fi,..., G C are the roots of p. By assumpfion we know fhat |^, | > M for each 
i = 1,...,d. This implies that 

/(z) =lnp(z) = lnp( 0 )-F^ln(^l-^^ (28) 

for any z such that |z| < |f|. Using the standard Taylor expansion for fhe principal 
branch of the logarithm, we obtain 

/ X \ n 1 xm 


13 












where Rn satisfies 


Rn 


00 


E 

m—n^l 


1 r 
mXr 


< 


1 \t\^ 

(« + l) “ (w + l)(l-‘?)' 


since q = \t\/M < 1. Noting that 


1 

m 


E 


i=l 


^ m 
i 


m! rfz’"'' ^ 


z=0 


(|2^ now follows by combining ll28ll with (l29ll and using the bound on | Kn |. 

Take n = C(lnd/e) (where C > (Inl/ij) ^is large enough so that 1/n < 1-q). 
Then the right-hand side of ll26ll is at most e. Write z = Then we have 

|e/(0-z| = |e/(d-z| < el/(f)-z| < gf similarly < e^. (This follows form the 

fact that for a complex number y = a + bi, we have \ey\ = e“ and \a\ < \y\.) Moreover, 
the angle between and is bounded by |51ne^ < \ lne^ < e. This shows 

To see ll24ll . take n = 0(ln 1/e) such that the right-hand side of ll26ll at most de. Then 
we have |e^| < and < |e^|e*, which is equivalent to ll24t . □ 


4.2 Approximating tensor network contractions 

In this section we will combine Theorem with Lemmas [TO] and [TT] to give a quasi¬ 
polynomial time algorithm for the computation of tensor network contractions for a 
large class of tensor networks. In particular. Theorem |2] will be proved here. 

Consider a graph G = {V,E) and h (z Sq. Suppose that h satisfies 




1 r(A+i) 

1.012(A(G)-kl) 


(30) 


for all e y and a e Let I = (F)z;ev G Sg be equal to the all-ones tensor for 

?ach V 

1^1 by 


each V e V, i.e., 1^ = 1 for all a e 


. Define a univariate pol 5 momial q of degree 
q(z):=p(G)(I + z(h-l)). (31) 


Note that q(l) = p(G) (h). So we are interested in computing the value of at z = 1. By 
Corollary!^ we have that for any |z| < 1.01, q{z) 7 ^ 0. So by Lemma ITT] we can compute 
a fast e-approximation to p{G){h) provided we can compute the mth derivative of q at 
z = 0 efficiently. This can be done in time for each m. Indeed, for any m, 


jm jm 

^'?(z)iz=o= E ^n(-f+z(^--f))(</’(^(^)))L=o 

<p-E^[k] veV 

, , 

= E E n(^-^)(</’(‘^(^)))^^'^'iz=o 

<l)-.E^[k]ucvveu 

= E E (32) 

UQV,\U\=m(p:E^[k] veU 


Let us denote by E (Lf) the edges in E that are incident with some vertex of U for U C V. 
Then (l32ll is equal to 

m! E 

UCV,lUI=m ip:E(U)^[k]veU 


and hence can be computed in time _ |y| 0 (m) gy Lemmas ITOl and [TTl we 

have the following result, which clearly implies Theorem |2l 
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Theorem 12. Let e > 0. Then for any graph G = (V,E) and any h ^ Sq such that h 
satisfies (l30l for all v E V and a. E zve can compute a multiplicative e-approximation 

to p{G){h) in time and an additive e\V\-approximation to In \ p{G){h) \ in time 

|y|0(ln(l/£))^ 

Remark 2. Let A e N and fix an edge-coloring model h such that h satisfies (l30l l for all 
d < A. Then by Corollary[6^ we have that \p{G){h)\ > for some constanf c > 0 
(assuming G is connected). So Theorem [121 provides a deterministic FPTAS to compute 
ln|p(G)(lz)| on graphs with maximum degree at most A. 

Remark 3. Theorem [12] of course also applies to tensor networks that strictly satisfy the 
conditions of Corollaries [6^ and b, but we leave the details to the reader. 

4.3 Approximating evaluations of graph polynomials 

In this section we will combine a result of Csikvari and Frenkel fT9l wifh Lemma [lU| 
and Lemma ITT] to give a quasi-pol 5 momial approximation scheme for evaluating a large 
class of graph pol 5 momials, including the Tutte polynomial. 

Let X ■ 0 C he a graph parameter. Define for a graph G = {V,E) the graph 
polynomial p;^.(G)(z) := Ei=i;\:/c(G)z^ where for k = 1,..., |y|, 

Xk{G):= ^ Ylx{G[Vi]) (33) 

y=ViU...ui4 !=i 
Vi C V nonempty 

y-nv,—0 if 

(here G[l^] is fhe subgraph of G induced by 1^). So x(G) is a polynomial of degree af 
most |y| with zero constant term and the coefficient of z is equal to x(G). If 7 ^ 0, 

fhen the degree of p^ is equal fo |y|. Csikvari and Frenkel |[T9l call graph pol 5 momials 
fhat arise this way graph polynomials of exponential type. (In fact, they give a different 
definition, but show that the definition above is equivalent to it.) 

A graph pol 5 momial of exponential type p(G)(z) = EEi cii{G)z^ is called of bounded 
exponential type if it is monic of degree n = |V(G)| for each G and if there exists a 
function R : IN —t [0, oo) such that for each f > 1 and v E V we have 

^ |ai(G[S])| <R(A)^-i 

scy(G) 

D€S,|S|=f 

Csikvari and Frenkel proved fhat these pol 5 momials have bounded roots on Q^: 

Theorem 13 (Csikvari and Frenkel 1191 ). Let p be a graph polynomial of bounded exponential 
type. For any A e N there exists a constant c > 0 (one may take c ~ 7.30319) such that for 
each G E the roots of p{G) are contained in the set {z e C : |z| < cR(A)}. 

Csikvari and Frenkel showed that many well known graph polynomials are of 
bormded exponential type such as the chromatic polynomial, the Tutte polynomial, 
the modified mafching pol 5 momial and the adjoint polynomial. Let us call a graph 
parameter x ■ G ^ G efficient if it can be computed in time We will show 

that we can approximate evaluations of graph polynomial of bounded exponential type 
using Lemmas [lO] and [TTJ provided x is efficient. 

Theorem 14. Let A e N and let G = {V,E) be a graph with maximum degree at most A. Let 
X '■ G ^ Cbe an efficient graph parameter with x{Ki) = 1 and suppose that all roots ofp^{G) 
have absolute value at most c,for some c > 0. Fix k > 1. Then for any e > 0 and x E C with 
\x\ > Kcwe can compute a multiplicative e-approximation to p^{G){x) in time 
and an additive e\V\-approximation to In \p^{G){x)\ in time |l/|0(lni/e)^ 
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Proof. Define the graph polynomial q by 

^(G)(z) :=zl'^lp;^(G)(l/z). (34) 

Note that q{G){z) = e|=[|X*:(G)zI^I“*^. Consequently, the roots of q, for G e G/^, are 
contained in the set A := {z e C | |z| > \}. Note that q(G)(0) = x\v\ = 1/ = 1- 

We will now show how to compute the derivatives of q. 

Let us write n := | V| and define for k = 0,1,..., n, Pj- to be the set of partitions of 
V into exactly k nonempty sets. Then 




z—0 


n k jm 

E E nx(GlVil)fs-'-^ 

feOIVi. 


For m n — k the contribution will be zero. So l l35ll is equal to 


(35) 


n—m 

E n^(G[^^;])- (36) 

{Vi,...,V„-m}SPr:-m 1=1 

No more than m sets in Vn-m can have size greater than or equal to 2. So to enumerate 
Vn-m, we first select a set of n — 2m vertices and then find all partitions into exactly m 
sets of the remaining 2m vertices. This gives a total of steps. Since x is efficient, 

the sum (l36l l can be computed in time 

Let now t = 1/x. Then |t| < Lemmas [TOl and ITT] then imply that we can compute 
numbers in time such that e~‘^ < \-^\ < P and such that the angle 

between and q{t) is at most e, and ^2 in time such that | In \q{t)\ — ^ 2 ] < £W- 

Then cci := is a multiplicative e-approximation to Pxi^) := ?2 + win|x| is 

an additive en-approximation to In | p;j-(x) |. This finishes the proof. □ 

Csikvari and Frenkel Remark 5.4] proved that the Tutte pol 5 momial is of expo¬ 
nential type (it is clearly multiplicative). It also follows from their characterization of 
graph polynomials of exponential type, cf. IfT^ Theorem 5.1] that the associated graph 
invariant x is efficient. Jackson, Procacci and Sokal l27l Theorem 1.2 and 1.3] showed 
that for every fixed uq G C there exists a constant C = C(A, uq) such that the roots of 
T(G)(q,Vo) for any graph G with maximum degree at most A are contained in a disk of 
radius C. This implies that Theorem |3] follows from Theorem [141 


5 Sparse graph limits and edge-coloring models 

In this section we will prove Theorem |4] To do this we will use the framework of 
Csikvari and Frenkel 1191 . 

Let A G N, and consider for any p > 0 the following collection of edge-coloring 
models 

£(A,q) :={h:N^^C: \h{ci) - 1\ < (l-p)E^±^ for all a G N^}. 

Recall that f*{A -t-1) > f*{l) ~ 0.71885; see Section |3] Now fix some small constant 
> 0 and fix h G £{A,S). Consider the following two graph polynomials: 

q{G){z) := p{G){I — z{I — h)) and q{G){z) := k~^^P'^^zPP'^^q{G){l/z). (37) 

Observe that p(G)(l) = p{G){h) = fcl^(‘^)lp(G)(l) and that q is monic. 
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Since for |z| < we have I + z{h — I) e £{A,0), it follows by Corollary | 6 ^ that 
p(G)(I + z(h — I)) 7 ^ 0 for all graphs of maximum degree at most A. This implies 
that q{G){z) 7 ^ 0 for all z e := {z e C : |z| < ^ 37 } and all G e Q/^, which in 
turns implies that for all G (z Q/^, all roots of q(G) are contained in the the compact set 
:= {z e C : |z| < 1 - J}. 

For a graph G we define }Iq to be the uniform distribution on the roots of if(G). 
Note that for all graphs G € Qa the measures }Iq are all supported on K^. We have the 
following result: 

Theorem 15. Let (Gn) be a locally convergent graph sequence. For any open set O c IR'^ and 
any continuous function / : x O — 7 IR that is harmonic on the interior of Kg for any fixed 

CO & n, and harmonic on Cl for any fixed z e X^, the sequence 



f{z,co)dyG„{z) 


converges locally uniformly in co G Cl to a harmonic function on Cl. 

We will prove Theorem [15] below. Let us first note that it implies Theorem |4] by 
proving the following more concrete result. 

Theorem 16. Let A e N and let S > 0. For any edge-coloring model h e £{A,5), the 
normalised partition function n{-){h) : Qa —7IR is estimable. 


Proof. Let (Gn) be a locally convergent graph sequence with Gn G Qa for each n. Take 
n = C \ X^ and let f{z,co) := In |a; — z| on X^ x O. Then taking 1 G O, we have by 
Theorem [15| that 

/ f{z,f)dyG„{z) (38) 

JKg 

converges. Now note that (l38ll is equal to 


I ln|l-z|dpG„(z) = 

JKs 


ln\q(G„)(l)\ 


\V{G„)\ 

ln|/c-l^(G")lp(G„,fe)| 

\V{Gn)\ 


\V{Gn)\ 


(39) 


Now we just need to observe that since (Gn) is locally convergent, we have that = 

converges. So (i39l l implies that the sequence {n(Gn){h)) is convergent, proving 
the theorem. □ 


Remark 4. Theorem [T^ is stated only for edge-coloring models that are close to the 
all-ones model. Of course it also applies to edge-coloring models that strictly satisfy 
the conditions of Corollaries |6^ and b. In particular, in view of Remark |TJ the result 
gives a qualitative version of a result of Borgs, Chayes, Kahn and Lovasz UTI on 'right 
convergence'; see also Il30l . We however leave the details to the reader. 

It now remains to prove Theorem ITSl 

Proof of Theorem \T5\ We will apply Theorem 4.6 from ||19|. Write q{G){z) = flo(G) -l- 
fli(G)z^ -|-... -|-fln(G)z” with n := |1/|. The only requirement in lT9l Theorem 4.6] that 
is left to check is that the coefficients «,,_, (■) can be expressed as a linear combination 
of the parameters ind(H, ■) for i = 0,...,m for each fixed m G IN. (In fact, Csikvari and 
Frenkel do not work with induced subgraphs, but with subgraphs, but these notions 
are equivalent by inclusion exclusion.) 

We will show that for any graph G = (V, E) G on n vertices the coefficients of 
z^,... ,z’” of q{G) can be expressed as a linear combination, over graphs H of at most 
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(A + l)m vertices, of the parameters ind(H, G). Since q(G) (z) = k a„_iz\ by 

dSZJ, this is clearly equivalent to the statement about the coefficients of q. 

By definition, q{G){z) is equal to 


E Yli^-iz{i-h)){(pi5{v))) = ^ E^'E n(^-^)(‘/’(^(^))- 


(p:E^[k]veV 


(p:E->-[k]i=0 UCV ueU 
\U\=i 


So we only need to show that we can express 

E L Ylii-h)imv)))= E E Uii-h){mv))) 


(40) 


(p:E — y\k] liCV u^il 
\U\^i 


(p:E — y\Jc]u^U. 

\U\—i 


, m. 


as a linear combination of the parameters ind(H, G) for each f = 0,. 

We need the concept of a fragment, which is a pair (H, i/^), where H is a graph and 
where ip is a map ip : V (H) —t {0,1,..., A}. We think of f>{u) as a number of half edges 
sticking out of u. For UCV we let G(U) be the fragment {G[U],f>) where f{u) is equal 
to the number of edges that connect u with V\U. Clearly, for each U of size i the second 
sum on the right in lIlOll only depends on the isomorphism class of the fragment G(U). 
(An isomorphism from a fragment {H,^p) to a fragment {H,ip') is an isomorphism a of 
the underlying graphs such that for each u G V{H), f>{u) = xpfa.{u)).) For a fragment 
F = {H,ip) let E(F) denote the set of edges of F including half edges. Then define for 
an edge-coloring model y. 


p{F){y)-= E n y(<?’(^(^)))- 

(p:E{F}^[k\v€V{H] 


Then we can rewrite ll40t as 


j- fc|E(G)|-|E(F)lp(p)(7_/j)ind*(F,G), 

F={H,<p) 

\V{H)\=i 


(41) 


(42) 


where the sum rims over fragments and where ind* (F, G) denotes the number of sets 
U of size \V{H)\ such that G{U) is isomorphic to F = 

So to show that (l42l l can be written as a linear combination of the parameters 
ind(H, G) for certain graphs H, it suffices to show that we can do this for ind* (F, G) for 
any fragment F. To see this, note first that we can write 

ind(H,G)= ind* {{H,4,),G). 

Next, say that xp < k for two maps ip,K : V{H) —t {0,1,..A} if xp{v) < k(v) for all 
V E V(H). This defines a partial order on := {xp : V(H) —)• {0,1,..., A}}. Let Z 

be the associated zeta matrix; that is, Z is the x matrix defined by. 




1 if Xp < K 

0 otherwise. 


Let a E be the vector defined by a^p = ind*{{H,xp),G). Then ind(H,G) = 

{Za)^g, where xpQ denotes the all-zero map. For any fixed xp E (Za)^ is equal 

to X]h' ind(H', G), where this sum runs over all graphs H' on | V(H)| -|- YZveV{H) fi'^) 
vertices containing H as an induced subgraph such that for each v E V{H) this vertex 
has exactly xp{v) extra neighbours in H'. It is well-known (by Mobius inversion) that Z 
is invertible. Hence a = Z~^{Za), which gives for each fragment {H,xp) an expression 
of ind* {[H, xp), G) as a linear combination of ind(H', G), where H' contains at most 
(A -|- l)m vertices. This finishes the proof. □ 
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6 Concluding remarks 

In this paper we built on ideas of Barvinok and Soberon J8l to find regions where the 
partition functions of an edge-coloring models does not evaluate to zero. We more¬ 
over used a technique of Barvinok |4j to transform these zero-free regions into quasi¬ 
polynomial time approximation algorithms for computing partition functions, as well 
as for evaluating graph polynomials. This work leads to some questions: 

• As is indicated by Example |2l our results are not optimal. Can one determine 
more precisely the region where the partition functions does not evaluate to zero 
on bounded degree graphs? Recent work of Barvinok fT) indicates that also zero- 
free regions that are not shaped like a ball are interesting. In [7] he used a strip 
argument to find a quasi polynomial time algorithm to approximate the perma¬ 
nent of a real matrix A for which <5 < Aj j < 1 for some fixed <5 > 0. 

• Can zero-free regions be used in any way to yield (randomized) polynomial time 
approximation algorithms? 
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